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Abstract We present a study on the dynamics of a system consisting of a pair of hardcore 
particles diffusing with different rates. We solved the drift-diffusion equation for this model in the 
case when one particle, labeled F, drifts and diffuses slowly towards the second particle, labeled M. 
The displacements of particle M exhibits a crossover from diffusion to drift at a characteristic time 
which depends on the rate constants. We show that the positional fluctuation of M exhibits an inter- 
mediate crossover regime of subdiffusion separating initial and asymptotic diffusive behavior; this 
is in agreement with the complete set of Master Equations that describe the stochastic evolution of 
the model. The intermediate crossover regime can be considerably large depending on the hopping 
probabilities of the two particles. This is in contrast to the known crossover from diffusive to subdif- 
fusive behavior of a tagged particle that is in the interior of a large single-file system on an unbound 
real line. We discuss our model with respect to the biological phenomena of membrane protrusions 
where polymerizing actin filaments (F) push the cell membrane (M). 

Keywords Brownian ratchet, driven dynamics, subdiffusion, non-equilibrium fluctuation 
1 Introduction 

Single-file diffusion of a system of hardcore particles is a well-studied process that provides a 
basic description of transport, for example, in fast ion transport through channels flO], in zeolites 
J^t] and in superionic or organic conductors 0), While their collective diffusion is like that of a set 
of independent particles, the diffusion of a tagged particle is known to be different - namely, its 
mean squared displacement (MSD) is proportional to \ft and its positions are Gaussian-distributed 
J5-12]. It is an exactly solved problem in the case when the system consists of identical particles 



having the same diffusion constant. On the other hand, the case when their motion is characterized 
by a set of different diffusion constants is not readily amenable to an exact mathematical analysis. 



However, Ambjornsson et ai, 11131,11411 have recently shown that the diffusion of a system consisting 
of just two hardcore particles with different diffusion constants on an unbounded (one dimensional) 
real line can be solved exactly; in particular, they have shown that the MSD of a tagged particle is 
proportional to i; they have also presented Monte Carlo evidence to show that it crosses over from 
diffusive (oc t) to subdiffusive (oc \/t) behavior only if it is in the interior of a large system. 

In this context, it is of interest to study the effect of constraining boundaries on the single-file diffu- 
sion of particles with different diffusion constants. Such boundary effects are relevant, for example, 
to understand the physical mechanism underlying the process of cell protrusion where polymerizing 
actin filaments push the cell membrane. The first physical description of the cell protrusion process 
was based on the Brownian Ratchet (BR) model II15L a. This is a one-dimensional two-particle 
representation for the filament-membrane system. In this model, the random diffusive motion of the 
cell membrane, represented as a Brownian obstacle, is rectified by the growing tip of a semi-rigid 
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rod whose other end is fixed at the origin. Using the stationary distribution for the gap between the 
tip of the flexible rod and the obstacle, they could explain the experimentally observed load-velocity 
curves reasonably well. 

In the two-particle model the membrane is represented by a single particle where the tension- 
induced correlations among the constituents of a two-dimensional flexible surface are neglected. 
In fact, a recent simulation study on two interacting random surfaces [17 [representing the cell 
membrane and the actin cortex has indicated the importance of correlation j l8|| . The representation 
of a long semi-rigid actin filament by a single particle is reasonable as long as nucleation of new 
filaments near the leading edge of the protruding membrane is negligible. In this case, a single actin 
filament polymerizes and depolymerizes at its both ends such that shrinking and growing of its length 
would never lead to dissolvation. This assumptionjs justified based on the observation that in vivo 
the ATP-mediated (de)polymerization processes 11911 can yield an effective 'treadmilling' 12011 of a 
filament which, on average, depolymerizes only at its 'minus' end and polymerizes only at its 'plus' 
end at the same rate such that its average length remains constant. Therefore, the theoretical ansatz to 
consider only the growing plus end, which is located near the protruding edge of the cell membrane, 
and to model this by a single particle is reasonable. Furthermore, in the two-particle model the 
assumption is made that the whole filament is prohibited to perform large-scale thermal motion 
(e.g. diffusion). In vivo this immobility is mostly achieved by strong adhesive bonds connecting the 
filament via the cell membrane to the underlying substrate on which the cell crawls. Therefore, in 
the present two-particle model the immobility of the filament is presumed, and the only mechanism 
by which a cell can be moved is by pushing the leading front of the cell membrane (particle M) by 
the polymerizing plus end of a filament (drifting particle F). 

We present an exact analysis of a two-particle system drift-diffusing on the positive real line in the 
case when particle F drifts and diffuses slowly towards particle M. The fluctuation in the position 
of the membrane particle M exhibits diffusive behavior at short, as well as at long times with an 
intermediate crossover-regime that is due to the onset of hardcore interaction between the particles. 
We show that this is in agreement with the results obtained from the set of Master equations that 
describe the stochastic evolution of the two-particle model. On the other hand, fluctuation in the 
position of the pushing 'filament' particle F is diffusive without any such crossover behavior. The 
mean displacement of the Brownian particle M ('membrane') crosses over from an initial diffusive 
behavior to an asymptotic drift behavior. For the sake of completeness we consider also the case of 
a repulsive boundary for F at the origin. 

In the next section, we describe the one dimensional drift-diffusive motion of a pair of hardcore 
particles in terms of a set of hopping rates on a regular lattice and present a discussion of this process 
in the continuum limit. In section III, we compare the results obtained by solving the two-particle 
drift-diffusion equation with those obtained by numerically integrating the set of master equations, 
given in the Appendix, that describes the hopping process on a lattice. 



2 Two hardcore particles on a Id lattice: 

Drift-Diffusion equation for the joint probability distribution 



Let a pair of particles, labeled F and M, be at positions np and tim respectively on a one di- 
mensional lattice; hardcore interaction ensures that np < n^j at all times, if that was the case at 
t = 0. Equivalently, the distance of separation between these particles, n = {nu — np), satisfies 
the inequality n > 1 at all times. A schematic illustration of the model is shown in FigQ] 

Consider a discrete time process. At every instant of time, particle M moves one step to the right 
or to the left with a priori probability pm or (1 — pm) respectively; on the other hand, particle F 



3 




FIG. 1: A schematic of a system of two hardcore particles in the presence of a reflecting boundary. In the 
discrete-time picture, the 'driven' particle, M, jumps right (or left) with a priori probability pm (or qM = 
1 — pm) at every instant of time. The 'driving' particle, F, jumps right (or left) with a priori probability pf (or 
qF — 1 — Pf) once in 1/wf instants of time on the average. Shaded region denotes the domain of two-particle 
drift-diffusion. 



may either move or stay put. Let wp be the probability that F will move at any given instant of 
time. Then, on the average, F moves just once during the period of time in which M has moved 
1/wp times. It is clear that wp also denotes the a priori probability that both F and M will move 
simultaneously at a given instant of time. Given that F moves, let it move one step to the right or 
to the left with a priori probability pp or qp = (1 — pp) respectively. Of course, with respect to 
the physico-biological situation, as described in the Introduction, the practically reasonable case is 
qF < pf, which implies rip > 0. 

So, when both the particles jump simultaneously in the same direction, the separation distance, n, 
does not change; when they jump simultaneously in opposite directions, n changes by ±2; on the 
other hand, when particle F does not jump, n changes by ±1. Let q^, qf, denote the probabilities 
per unit time that the change in separation distance 5n = 0, ±1, ±2 respectively. They are given by 
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= wfPfPm 




= w F {\ -Pf){1 - Pm) 


it 


= (1 - Wp)p M 


<h 


= (1 - w F )(l - p M ) 




= w F (l -p F )p M 


<h 


= w f pf{1 - Pm) 



and they all add up to unity. 

The physically motivated constraint that the particles can only move on the positive real line, 
7i m > np > 1 (on a lattice), as well as the probability, wp, for them to move simultaneously imply 
that we have different sets of Master Equations for the joint probability, P(np, um', t) corresponding 
to the cases n = 1, 2, > 3 respectively (see, Appendix). 

In the continuum description, the joint probability for the positions of the hardcore particles, xf 
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and xm, satisfies the following drift-diffusion equation, 

dP(x F ,x M ]t) ( d 2 P{x F ,x M ;t) dP(x F ,x M ]t) 



dt \ F dxp ^ F dxp 

d 2 P(x F ,x M ;t) dP(x F ,x M ;t)\ 
D m 5-3 Mm t— I (2) 

ox M oxm J 

which is, in fact, the continuum version of the unrestricted Master Equation, Eq( l36l l. Hardcore 
interaction between the particles constrained to be on the positive real line implies that xm > x F > 
at all times. In our model, particle M is purely diffusive whereas particle F can drift as well. 
Therefore, as shown in Appendix (f), the corresponding diffusion and drift coefficients are given by, 

D F = -ji.lt + <h + <l0 + %) 

±w F (l-w F ) + 2w F p F (l-p F yj ~^w F (w F small) 



y>F = -?2 + <?2 + % - % = w f(2PF - 1) 



d m = -j[qt + i2 +<ii + <ii +<io+%) = \ 



(3) 



Hm = q£ ~ ? 2 +Qi -Qi + Qo ~ % = 
The initial condition for this problem could be specified by the joint probability distribution, 

P(x F ,x M ;t = 0) = S(x F - x° F )6{x M - x° M ) (4) 
where S(x) is the Dirac delta function. The first boundary condition, 

dP(x F ,x M ;t) dP(x F ,x M ;t)\ 

D ™ °F =l(»M-VF)P(x F ,X M ;t)} Xf=XM (5) 



8xm dx F 



expresses the fact that the particles cannot pass each other (x F < xm). The second boundary 
condition, 

dP(x F ,x M \t) \ 

Up fi F P(xp,x M ;t) =0 (6) 

dx F J XF=0 

ensures that there is no current across the boundary at x F = 0. These two boundaries define a wedge 
as the domain for the joint distribution P(x F , xm t) ■ The question is whether a separable solution to 
Eq.© can be found. We first transform the variables, {x F , %}, into a pair of 'collective' variables, 
say {x, r}. 

Inter-particle separation, r = xm — x F , could be fixed as one of the new variables. For the 
other variable, we set x = cmXm + c F x F and try to fix the dimensionless constants, cm and c F , 
by requiring that the transformed drift-diffusion equation also has the same form as Eq.©. This 
requirement leads to the condition that the coefficient of d 2 P(x, r; t) / dxdr vanishes: 

D m cm - D F c F = 0, (7) 

This implies that 
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where c is an arbitrary constant. We set c = 1/2 so that x will be the position of the center- 
of-mass of the system when the particles are of the same mass and have the same value for their 
diffusion constants {Dm = Dp). In general, x is not the center-of-mass coordinate because the 
equality, cm + Cf = 1, will be satisfied only when the particles have the same mass and also when 
Dm — Dp- Nevertheless, with the choice c = 1/2, we have the transformation, 



1 / D M D f 

X = z fe^fe 1 " 



r = xm - x F 



(9) 



under which Eq.© becomes 

dP{x,r;t) _ ( d 2 P(x,r;t) dP(x,r;t) 
dt \ X dx 2 x dx 

with the initial condition, 



D, 



d 2 P(x 1 r;t) dP(x,r;t) 



dr 2 



dr 



(10) 



P{x,r;t = 0) = 



(D F + D 



where 



\2y/ D F D M 



5(x - x°) 5{r - r°) 



(11) 
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(12) 



° F x° 

d^ Xm 



The boundary condition, Eq.©, that ensures hardcore repulsion between the particles now trans- 
forms into the following condition, 



D r 



dP 
dr 



P 







(13) 



r=0 



which is the reflecting boundary condition at r = (i.e., xp = xm) for the above drift-diffusion 
equation for P(x,r;t). The other boundary condition, Eq.©, transforms into 



1 , OP 



D 



dP 
dr 



VfP 







(14) 



It is clear that the wedge in {xf, XM}-space has transformed into a wedge in the {x, r}-space. 
Moreover, the above boundary condition implies that separable solution to Eq.dTob is not possible in 
the transformed space. 



A. No reflecting boundary at xp = 0: 



6 




FIG. 2: Shaded area is the domain of the joint probability distribution, P(xf,Xm', t), with only one reflecting 
boundary corresponding to the hard-core repulsion between the particles. 



One way to ensure separable solution to Eq.dTOt is to ignore the boundary condition, Eq.©, and 
hence Eq.dPfli. This amounts to assuming that the growing filament does not degrade to a monomer 
in any arbitrary time span of interest; so, particle F moves freely as long as it is to the left of particle 
M. In this case, we may now try a product solution to Eq.dTol) of the form, 

p ^ t) = {lm) p{x '' t)p( '' t) <15) 

so that P(x; t) and P(r; t) satisfy the equations: 



dP(x;t) d 2 P(x;t) _ ^ dP(x;t) 

dt x dx 2 x dx 



dP(r;t) d 2 P(r;t) dP{r;t) 

dt T dr 2 1 dr 



(16) 



The first one is a free-boundary equation with the initial condition, P(x;t = 0) = S(x — x°), 
whereas the second one, with the initial condition P(r;t = 0) = 8{r — r°), is subject to the 
boundary condition given by Eq.(fT3l). Solutions can be written down immediately jilll : 



P(x;t) = g 



(17) 



P{r , t ) = g( r - r l±_^ ) + e^°/»r g ( r + r°+Hrt 



/V ^/c, f°° f y + r-n r t 



Dr J y=r a \ y/AD r t 

where <?(£/ct) is the normalized gaussian function given by 



dy 



g{ii/a) = -4= e- e/ ° 2 (18) 
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The distribution P(r; t) will become asymptotically stationary when fi r > 0; this implies the condi- 
tion wf(2pf — 1) > (2j>m ^ 1) which is to say that the effective rightward drift of particle F should 
be more than that of particle M. It is clear from the definition, Eq.((9]l, that stationary value for (r) 
does not imply stationary values for (xm) and (xf). 

B. Reflecting boundary at x = 0: 




FIG. 3: Shaded area is the domain of the joint probability distribution, P(xF,XM',t), with one reflecting 
boundary at x — and the other corresponding to the hard-core repulsion between the particles. 



In the Brownian Ratchet model 1 15. 1611 a reflecting boundary at the origin is introduced in order 
to provide a load which prohibits backflow of the actin filament. In this model the (de)polymerizing 
filament is represented by a rod where the plus end can (de)polymerize and the minus end of the 
filament is fixed at the origin. Since particle F represents the plus end of the rod, which has to be of 
nonzero length, we have a reflecting boundary for particle F at the origin. 

This reflecting boundary at xp — transforms into the boundary condition Eq.(fT4b. which rules 
out a separable form for P(x, r; t). On the other hand, a separable solution to Eq.([T0l) can still be 
obtained if we arbitrarily impose the following boundary condition at x = 0: 



D 



dP(x;t) 
dx 



H x P(x;t) 







(19) 



2=0 



From the definition, Eq.©, we see that a; = implies .tj? = —{Dp / Dm)xm\ hence the domain for 
the joint distribution, P{xp, xm', t), is the wedge schematically shown in Figj3] It is unphysical if 
we insist that particle F represents the tip of a polymer rod; yet, for small wf (by definition, Eq.fO, 
{Dp /Dm) ~ wf), we may expect to have an approximate solution to the model. 

Assuming se para ble form for P(x, r; t) (Eq.(fT5ll), the above boundary condition at x = leads 
to the solution, 12111 



P{x: 



= fx-x°-y x t\ „ /> ('■ ■<■" /'■■') 



(20) 



The distributions, P(x;t), will become asymptotically stationary when fi x < 0. From the defini- 
tions, Eqs.([3lT2]i, we see that fi x = wf{pm+Pf — 1) and so /i^ < implies the condition pp < qM- 
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This will be realized when the leftward drift of particle M is more than the rightward drift of particle 
F. Interestingly, stationary value for (x) implies, by definition Eq.©, stationary values for (xm) and 
(x F ) also, which in turn implies stationary value for (r). Hence, the condition fj, x < is enough to 
ensure that the joint probability distribution, P(x, r; t), becomes stationary. 

C. Tagged particle distributions: 

Transforming back to the variables, x F and xm, with appropriate Jacobian prefactor, the above 
solutions for P(x; t) and P(r; t) lead to the joint distribution P(x F , XM',t) and hence to the tagged 
particle distributions: 



/>oo 

P(x F ;t) = / dx M P(x F ,x M ;t) 

P(xM',t) = / dx F P(xF,XM;t) 
J —oo 



(21) 



The upper limit x^j suggests that x F may be arbitrarily close but never equal to xm\ similarly, the 
lower limit suggests that xm may be arbitrarily close but never equal to x F . The Mean Squared 
Displacement (MSD) of a tagged particle, a 2 M F , may then be obtained from its corresponding dis- 
tribution. 



3 Results and Discussions 



A. Inter-particle distance, r: 

Stationarity for the distribution P(r: t) in Eq.dTTl) is ensured by the condition /i r > which, from 
Eqs.( l3ll2"b . implies w F (2p F — 1) > (2j>m — 1). In this case, the asymptotic stationary form of P(r) 
is given by 

P(r) = P(r; t -> oo) = ^ e -^ r/Dr (22) 
which leads to the asymptotically stationary values for the moments: 

fJ-r 

(23) 

(r 2 )^oo = 

VMr 

Initial time-dependence of these moments and their approach to these stationary values can be ob- 
tained from their full forms: 

(r) = ^erfc(r°-) + i(r -M)erfc(-r°-)- 

D r ,,^ /ri . / n_Lv / D r t 



^ e^-°/^ erfc(r o +) + !dll e -(r°-f (24) 

2fl r V 7T 
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( r 2) = (ElY [(erfc(r -)- e ^°/^erfc(r 0+ ) 



2D r 

Mr 



^ e -(-°-) 2 - -(r° + / M)erfc(r°+) 

7T 2 



i [(M) 2 + 2(D r » r r )t + (r ) 2 ] erfc(-r°-) + (r° - ^/^V^" (25) 



where we have the definitions, 



0- = ; (26) 

2 r _ 2 

erfc(x) = 1 = / e y dy 

V 71 " ./o 

The 'early' time behavior (t < r°/fj, r ) of the moments can be guessed by recognizing that r 0± ~ 
r° /y/4D r t is large and therefore, in the limiting case fi r — > for example, the average separation 
(r) increases proportional to \/t\ larger the value of /i r , shorter will be the growth regime (t) oc \/t 
because it will attain its stationary value D r //j, r faster. 

It will be interesting to compare these results with those obtained by numerically integrating 
the Master equations, Eqs. [38l l. given in the Appendix (e). We have presented in FigHJa) the 
distributions of separation-distance obtained by numerically integrating these equations for 133 
and 3398 time steps respectively; the distribution is clearly exponential at t = 3398 - namely, 
P(r) ~ 0.1e _0 096r suggesting that the average value (r) ~ 10.42. For the parameters used in the 
computation, the stationary value of (r) should be D r //j, r = (1 + wf)/2u>f = 10.5 according to 
Eq.d23b. Even the normalizing constant (~ 0.1) obtained from the Master equations is close to the 
expected value of fi r /D r ~ 0.095. 

In fact, the diffusion coefficient, D n , obtained from Eq.d3~8b in the continuum limit for the case 
?? > 3 is equal to D r (Eq.([T2b forjjjvf = 1/2 and wp small; but, /i rl = /j, r whatever be the value of 
Pm- 

The small difference in the approach of (r) to its stationary value, observed in Figj4fb) for the 
case p^j = 1/2, implies that the distribution obtained by integrating the Master Equations is not 
completely stationary even after 3398 time steps for the lattice size chosen; moreover, finite size 
effects could make slight differences because of the need to have boundary equations for n = L and 
n = (L — 1) similar to what we have for n = 1 and n = 2. 

B. Two-particle variable, x: 

In the case when there is no reflecting boundary at xp — 0, the variable x is unbounded and hence 
its distribution P(x; t) is a gaussian (Eq.(TP7li) with the following moments: 

(x) = x + flxt 

(27) 

(x 2 ) = 2D x t + (x + fi x t) 2 



fi r t 



(J, r t 



On the other hand, when the variable x is constrained to be positive, its distribution P(x; t) is given 
by Eq. d20l) , The corresponding moments can be obtained from Eqs. (1241 and ( T2~5l l by making the 
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following replacements: 



r x ; r° o x° 



/i r <-> -/x x ; D r ^ D x (28) 



The moments have the following asymptotic behavior {pL x > 0): 

(*>t-oo = M+f^ + ^e"^ /^ 
V Ma 

(29) 

(a; 2 ) t ^oo = (nxt) 2 + 2{D X + n x x° + D x e-»* x ° /D *)t + const. 

C. Mean Squared Displacement of a tagged particle: 

It is tempting to see how the mean-squared fluctuation in the position of a tagged particle, say M, 
would be related to those corresponding to the variables, r and x, given the fact that the transforma- 
tion (xm-i x f) < — > ( r i x) is one-to-one. Quite surprisingly, the following ansatz, 

°M = ( X m) - ( X M) = 77; , n so K + TTT a r (30) 



(D M + D F y \ x 4 D F 

with a x r calculated using the distributions, P(x;t) and P(r;t) given by Eq.(|20"l) and Eq.([T7]i re- 
spectively, leads to a striking agreement with those obtained from the Master equations. 

In Fig[2a), we have presented <tm obtained by using the above ansatz for the case when we have 
a reflecting boundary at x = (continuous line A) as well as from the Master Equations (open 
circles). The parameters are wp = 0.05;pp = 1;pm = 1/2 and the agreement is quite good. In the 
same figure, continuous line B represents <jm obtained using the above ansatz but for the case when 
there is no reflecting boundary at xf = 0. Its long-time deviation from both the case corresponding 
to a reflecting boundary at x = (continuous line A) and the Master eqns.(open circles) indicates 
that it does not represent the basic phenomenology of the Master eqns. 

On the other hand, a similar ansatz for the average positions 




M = 2 ^^[{x) + -J^{r)\ (31) 



Dm + Df 



does not lead to such a good agreement with what we obtain from the Master Equations at early 
times. However, it leads to qualitatively similar asymptotic behaviors - namely, the drift velocities 
are proportional to each other, the proportionality constant being Wp -dependent. For example, in 
(FigHb)), the average velocity of particle M obtained from the Master Equations (open circles) 
is presented along with that obtained using the ansatz Eq.(|31~b (line M). It is quite clear that their 
asymptotic drift velocities are just proportional to each other even though their early time behaviors 
are different. Line F in the figure is the velocity of particle F obtained by using the above ansatz. 
While it takes some time, say t m that depends on wp and pp , for M to start drifting with constant 
velocity, F starts drifting almost from the beginning. 

In Figjfja), we have presented <jm obtained by using the above ansatz for the case when we 
have a reflecting boundary at x = but with wp = 0.1 (continuous line); the agreement with the 
Master Equations' data (open circles) is reasonably good. On the other hand, in the case Wp = 0.2 
presented in Fig(6|b), there is no agreement between the analytically computed data (continuous 
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line) and those obtained from the Master equations(open circles). This demonstrates, as mentioned 
earlier, that results obtained with an unphysical reflecting boundary at x = may agree with those 
obtained with the physically meaningful reflecting boundary at xf = for small values of Wf 
(namely, wp < 0.1 withp^ = 1 andpM = 1/2). 

In Figj7ta), we have shown how the crossover-behavior of <jm contrasts with the monotonously 
increasing behavior of a p. It is only when particle F has no effective drift towards particle M that 
such a crossover is non-existent (see Fig|7jb)); larger the effective drift, sooner is the crossover seen. 
It remains to be answered whether it is due to the reflecting boundary at the origin. 

We have shown in Fig[8]the fluctuation in the position of M for the cases, with as well as without 
a reflecting boundary at x — 0. The parameters used are mentioned in the caption. It is clear that 
the crossover becomes more pronounced and shifted to later times in the presence of a reflecting 
boundary than in its absence. In fact, the crossover is due to the smaller diffusivity of particle F than 
that of particle M (Dp < Dm)- The effective drift of particle F leads to the saturation of the average 
separation-distance. 



4 Summary and Conclusions 



It is appropriate, at this juncture, to compare this model with the standard BR model II15LI16H . The 
ratchet mechanism, in the BR model, is due to a monomer squeezing itself in the gap between the 
filament-tip and the barrier particle (i.e., the membrane). The resulting growth of the filament-rod is 
impeded by the inward (forced) drift of the barrier particle. This leads to a steady state in which the 
average gap betwen the filament tip and the barrier remains constant. The average ratchet velocity, 
in the steady state, is proportional to the net polymerization rate of the filament and is given by 



vbr = o [a exp 



P (32) 



where a and (3 are the polymerization and depolymerization rates respectively; D is the diffusion 
constant of the barrier particle, which drifts towards the filament-tip under the influence of the di- 
mensionless force, //fc^T; and, 5 is the size of the intercalating monomer. 

In our model, the BR-parameters (S, a, (3, D and //fcgT) are all lumped into the parameters Wf 
and pp so that the drift coefficient fj,p, defined in Eq.(f3]l, is equivalent to the ratchet velocity, vbr- 

Hp = v B r (33) 

Stationarity for the gap distribution is ensured by the condition fip > fi^j, which is equivalent to 
the condition wp > in the case when hm = and pp = 1 (see also, FigHfb)). In other words, 
the effect of the load-force and the consequent inward drift of the barrier particle in the BR model 
is mimicked in our model by the parameter wp > 0, which also tunes the steady state value of the 
average gap, (r) (Eq.(|23l. So, our model can be thought of as a variant of the BR model 1 15, 3 
that addresses the positional fluctuations of the particles M and F as well. 

In summary, we have discussed a two-particle model of cell protrusion namely, a system of two 
hardcore particles diffusing with different rates. The Brownian particle, labeled M, experiences a 
random hardcore 'push' by another particle, labeled F; M may be referred to as the 'driven' particle. 

We have solved the equations in the case when the drift-diffusion of particle F is small and ob- 
tained exact expressions for the tagged particle moments. Physically, this corresponds to the situa- 
tion when the effective drift of particle F towards particle M is small. Since particle M represents a 
'membrane', fluctuations in its position is of interest when it is being 'driven' by another particle F. 
Computed from the exact solutions (and using the ansatz Eq.(l30il). our theory exhibits a crossover 
from an initial diffusive behavior to an asymptotic diffusive behavior. The existence and duration 
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of the crossover regime depends on the diffusivities of the particles. It compares very well with 
the fluctuation data obtained by numerically integrating a complete set of Master equations that de- 
scribe the stochastic dynamics of this model. Asymptotic diffusion of a tagged particle is, of course, 
expected for a finite system. 

Appendix - Master Equations for the two-particle system 

Let P(n F , tim] t) denote the probability that, at a given instant of time t, the particles F and M 
are at positions n F and um respectively. The distance between them, n = um — n F , satisfies the 
inequality n > 1 at all times; the jump probabilities per unit time, {q}, are defined in section 2, 
Eq.©. 

(a) Case, n = 1: 



dP(l,2;t) 
dt 



qiP(l, 3; t) + <fo P(2, 3; t) - (q+ + % + )P(l, 2; t) 



(34) 



dP(n F , n F + 1; t) 
dt 



q 2 P(n F - 1, n F + 2; t) + q 1 P{n Fl n F + 2; t) + 

q$ P(n F + 1, n F + 2; t) + qQP(n F — 1, n F ; t) — 
(q+ + q+ + q+ + q^)P(n F ,n F + l;t); (n F > 2) 



(b) Case, n = 2: 



<9P(l,3;i) 
dt 



q^P(l, 4; t) + q+P(l, 2; t) + g "P(2, 4; t) 
fat +Qi +<Z + )P(1,3;<) 



(35) 



dP{n Fl n F + 2;t) 
dt 



q 2 P(n F — 1, n F + 3; t) + q x P(n F , n F + 3; t) + q^P(n F , n F + 1; t) + 
qQP(n F + 1, 7iF + 3; i) + q$P{n F — l,n F + 1; t) — 
(32" + 9i" + 9i" + So + %)P{n Fl n F + 2; i); (n F > 2) 



(c) Case, n > 3: 



dP(l,l + n;t) 
dt 



^ P(2, n; t) + q[P{l, 2 + n;t) + 
q+P(l,n;t) + qoP(2,2 + n;t) - 
(V2 + 4 + ?r + Qo ) p (!= 1 + »; *) 



(36) 



dP(n F , n F + n; t) 



q 2 P(n F + 1, iif + n — 1; t) + q 2 P(n F — 1, n F + n + 1; t) + 
P(n F , n F + n — 1; t) + q x P{n Fl n F + n + 1; t) + 



q$ P(n F — l,n F + n — l;t) + q P(n F + l,n F + n + l;t) — 
[it +<h +1i +<h +Qo + qo)P(n F ,n F +n;t); (n F > 2) 
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These equations could be solved for P(n F , tim', t), subject to the initial condition P(n F , tim) = 
S np n o 5 nu n o at t = 0. Since n F > 1, we also should have P(n F = 0, um', t) = 0. 

(d) Tagged particle distributions: 

The tagged particle distribution functions, P F (n F ; t) and Pm{tim', t), can then be obtained from 
P(n F , tim; t) using the following definitions: 

oo 

P F (n F ;t) = P(n F ,n M ;t) 

(37) 

TIM— 1 



PM(n M ;t) = 2_j P(n F ,n M ;t) 



np — l 

(e) Master Equations for the separation-distance, n, between the particles: 

The probability distribution function, P(n; t), for the separation distance, n, can be obtained from 
the above set of equations for the cases n = 1, 2, 3, ■ • • by summing over n F and also ignoring terms 
that violate the constraint n F > 1. 



dP(l;t) 

at 

dP(2;t) 
dt 

dP(n; t) 

at 



q^P(3; t) + q^P(2: t) - (q+ + g+)P(l; t) 
q^P(A; t) + P(3; t) + q+P(l;t) - 

{ql +lt + qi)P(2;t) (38) 

q^P(n + 2; t) + q+P(n -2;t) + 

q^P(n + 1; t) + qfP{n - 1; t) - 

(qt + ?a + qt + ?f *); ( n > 3 ) 

When summed up, these equations lead to the expected conservation of the total probability namely, 
S^Li Pn(n; t) = c where cis a constant. 

(f) Drift and Diffusion coefficients for the case, n > 3 and n F > 2 (second of Eq.( f36b ): 

Using the continuum variables, x F and xm, for the positions of particles F and M, we can rewrite 
the second of Eq.(l36b in discrete time as a jump equation, 

P(x F ,x M ;N + 1) = q£P(xF+l,XM-l;N) + q2P(x F -l,XM + l;N) + 

q+P(x F ,x M -l;N) + q^P(x F ,x M + l;N)+ (39) 
q£P(x F - I, x M -l;N) + qoP(x F + I, x M + I; N) 

where I is the infinitesimal displacement of the particle. Taylor-expanding the RHS of the above 
equation, and ignoring I 3 and higher powers, we have 

dP 1 d 2 P 

P(x F ,x M \N+l) « P(x F ,x M ;N) - [w F (2 PF ~ 1)1} j^— + - w F l 2 ^-j- 



dx F 2 dx ± 
BP 1 B 2 P 
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with the initial condition, 

p(x F ,x M ;0) = S(x F )S(x M - I) (41) 

For simplicity, we can take x F and x m to represent the displacement of the particles from their 
respective initial positions. Then the initial condition becomes 

p(x F ,x M ; 0) = S(x F )S(x M ) (42) 

In the absence of the reflecting boundary at x F = and the hard-core constraint, the variables x F 
and xm are unbounded, and we define the Fourier Transform, 

/oo />oo 
dx F e ikFXF I dx M e ikMXM P{x Fl x M ]N) (43) 
-oo J — oo 

From the initial condition, Eq.(l42b. it is clear that P(k F , kM', 0) = 1. Fourier- Transforming Eq.(l40t. 
we get 



1 



P(k F ,k M ;N + 1) w {1 + i[w F (2p F - l)l]k F - -w F l 2 k 



■2 1,2 

F 



+ i[(2 PM - l)l]k M - \l 2 k 2 M }P(k F , k M ;N) (44) 

which becomes 

iY 



P(k F , k M ; AT + 1) « ( 1 + i[w F (2p F - l)l]k F - ^w F l 2 k F + i[(2p M - l)l]k M - -l 2 k 2 M 



^w F l 2 k 2 F + i[{2p M - l)l]k M - 

'(45) 

when iterated with respect to the discrete time variable N and subject to the initial condition, Eq.d42l. 
Written in the equivalent form, we have 

P(k F , k M ; N + 1) w cxp(iVlog[l + a F + a M }) (46) 

where 

a F = i[w F (2p F — l)l]k F — —w F l 2 kp (47) 
a M = i[(2p M - l)l]k M ~ -^ 2 k 2 M (48) 
Using the standard expansion for the logarithmic function, 

oo n 

log(l +x) = V (-l) n+1 — ; (x small) (49) 



we expand P(k F ,kM', N + 1) given by Eq.d46l): 

P(k F ,k M ;N + 1) w exp(JV{a F + ajw + 

ii4(2 PJ r - i) 2 z 2 4 + ^(2pm - i) 2 / 2 *m + 

w F (2p F - \){2p M - l)l 2 k F k M }) (50) 
where we have ignored kp, k M and higher powers. With the following notational simplifications, 

p, F = w F (2p F — 1)1 (51) 
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HM = (2 PM - 1)1 (52) 



( ^w F (l - w F ) + 2w F p F (l - PF ) ) I 2 



D M = 2 PM (1 - p M )l 2 (54) 

we have, 

P(fcp, k M ; N + 1) ps exp(N{i^L F kF - D F k 2 F + i^m^m - D M k 2 M + MFMAffe.F&Af }) (55) 



which upon Fourier inversion gives 



^ w,tii °Hw)'' K1 <56) 



where 



= a; M - ^mN 

^,2 ,,2 

'Af = U M 



M 



Dt.r = L>M — 

AD F 



and is defined in Eq.dTSTl. It is clear that P(xf, xm', N + 1) will have a separable form only 
if hf = or = or both. In that case, the /x's and the D's defined in Eq. (l51H54l i will be 
corresponding drift and diffusion coefficients. 
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FIG. 4: (a) Probability distributions for the separation distance, r. Lattice size, L — 1024; at every time step, 
particle M moves left or right with equal a prior probability 1/2 whereas particle F moves always to the right 
(Pf = 1) with probability wf = 0.05. Open circles represent data dumped after t = 133 time steps (ti) and 
also after t = 3398 time steps (t 2 ). Straight line fit to the data (t 2 ) indicates that P(r) ~ 0. l e -° 096r after 3398 
time steps, (b) Average separation distance - open circles (Master equations); filled stars (analytic expression). 
The initial slope before saturation is (r) ~ \ft. 
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FIG. 5: Fluctuation for the tagged particle M, ayi = {{x 2 M ) — (in;) 2 ) 1 ' 2 . Master equations (open circles); 
ansatz Eq.J30t (continuous line). Parameters are the same as in Fig|4]- namely, (a) wf = 0.05; pf = 1; Pm = 
1/2. (b) Average velocities of the particles, M and F, obtained by using the ansatz, Eq.lBlt (lines M and F); 
they approach the same value asymptotically, implying that the particles drift together as a single system. Open 
circles represent the average velocity of particle M obtained from the Master Equations. It is clear that the 
asymptotic drift velocity is just proportional to that given by the ansatz. Parameters are the same as in Fig|4] 
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FIG. 6: Fluctuation for the tagged particle M, gm = {{ x \i) — (scm) 2 ) . Master equations (open circles); 
ansatz Eq.J30t (continuous line). Parameters: (a) wf = 0.1; pf = 1;pa/ = 1/2; (b) wf = 0.2; pf = 
1-Pm = 1/2. 



20 



10 1 - 



1 1 1 — I — I — I I I | 1 1 1 — I — r 

Reflecting boundary at x = 0; 


i i i | 1 1 1 — i — i i i i | r 

■ 


. (w =0.05; p = 1; p M = 1/2) 






















B : 'membrane' particle; 




B p : 'filament' particle. 
■ ■ ■ i 



10 1 



10' 



10 J 



10 4 



(a) 




(b) 



FIG. 7: (a) Computed fluctuation data (reflecting boundary at x = 0), Bm and Bf, for tagged particles M and 
F respectively. It is clear that ctm ~ Vt for Bf, whereas it shows a crossover from an initial diffusive behavior 
to an asymptotic (t ~ 5000) diffusive behavior. Parameters are u>f = 0.05;p/ = 1;pm = 1/2. (b) The ratio 
ffM/t 1 ^ 2 for different values of pf (= 0.5, 0.6, 0.7, 0.8, 0.9, 1.0; top to bottom) with wf having a fixed value 
0.05. For pf = 1/2 (i.e., no effective drift towards particle M), there is no crossover. 
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FIG. 8: Fluctuation in the position of M for the cases with (B) and without (A) a reflecting boundary at x = 0. 
The parameters are, wf = 0.05; pf = 0.9; pm = 0.5. The crossover is more pronounced for the case B than 
for the case A. The diffusivity of particle F, (Df = Wf/2), is smaller than that of M (Dm = 1/2). 



